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Abstract 

We investigate monoenergetic gamma-ray signatures from annihilations of dark matter com- 
prised of Z\ the first Kaluza-Klein excitation of the Z boson, in a non-minimal Universal Extra 
Dimensions model. The self-interactions of the non-Abelian Z l gauge boson give rise to a large 
number of contributing Feynman diagrams that do not exist for annihilations of the Abelian gauge 
boson B 1 , which is the standard Kaluza-Klein dark matter candidate. We find that the annihi- 
lation rate is indeed considerably larger for the Z 1 than for the B 1 . Even though relic density 
calculations indicate that the mass of the Z 1 should be larger than the mass of the B 1 , the pre- 
dicted monoenergetic gamma fluxes are of the same order of magnitude. We compare our results to 
existing experimental limits, as well as to future sensitivities, for image air Cherenkov telescopes, 
and we find that the limits are reached already with a moderately large boost factor. The realistic 
prospects for detection depend on the experimental energy resolution. 
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I. INTRODUCTION 



There is compelling observational evidence for a dark component of the matter content 
in the Universe, and yet, the particle nature of this dark matter (DM) remains unknown. 
Since there is no suitable particle DM candidate in the Standard Model (SM) of particle 
physics, the need for DM points towards physics beyond the SM. The most popular class of 
particle DM is weakly interacting massive particles (WIMPs), i.e., particles that are weakly 
interacting only and have masses in the GeV to TeV range. The standard example of such 
a DM candidate is the neutralino in supersymmetric extensions of the SM 

WIMP DM candidates also appear naturally in certain extra-dimensional extensions of 
the SM, i.e., so-called Kaluza-Klein (KK) theories. In particular, in the Universal Extra 
Dimensions (UED) model, which was proposed by Appelquist, Cheng, and Dobrescu in 
2001 jsj, the lightest Kaluza-Klein particle (LKP) is stable, due to conservation of KK 



parity. In the minimal version o 
B l of the U(1)y gauge boson B 



the UED model, the LKP is the first KK excitation 
3|, which is non-baryonic and electrically neutral, and 
hence, a possible DM candidate. This so-called Kaluza-Klein dark matter (KKDM) was 
first suggested in Refs. j^, 5|. According to detailed calculations of the relic abundance of 
B l DM, including the effects of coannihilations, the B 1 should have a mass in the range 
from 500 GeV to 1600 GeV in order to account for the observed DM . Taking; resonance 
effects from second KK level particles into account, it has been shown that the preferred 
mass range could increase to above 1 TeV 

The phenomenology of B 1 KKDM has been quite thoroughly investigated in the litera- 



14 



20]. 



ture [8-[l8|]. However, as shown in Ref. 19j], non-minimal UED models also support other 
LKPs, which could be DM candidates. Such scenarios give rise to different phenomenological 
implications, which have only been partially studied in the literature fiol . 

In this paper, we investigate the annihilation of KKDM into pairs of high-energy photons. 
This process results in monoenergetic gamma-rays with energy ~ itilkp, which is con- 
sidered to be a smoking-gun signature for DM. We study a non-minimal UED model, where 
the main difference to the minimal UED model is that boundary localized terms change the 
mass spectrum so that the LKP is the Z 1 . The flux of monoenergetic gamma-rays within the 
minimal UED model was calculated in Ref. 2l|. Furthermore, in Ref. 17j, the calculations 



were repeated, with similar numerical results, and in addition, the contributions to the mo- 



noenergetic flux from the final states 7Z and were calculated. These contributions were 
concluded to give only a small correction to the total flux. As a check of our calculations, 
we have successfully reproduced the analytical results in Ref. [21]. 



Gamma-rays from DM anni 
borne experiment Fermi-LAT [ 



filation are currently being searched for by the satellite- 
as well as by a number of ground-based image air 



Cherenkov technique (IACT) telescopes, such as H.E.S.S. |23|, MAGIC |2J], VERITAS |25|, 



and CANGAROO-III 



261 ] . Fermi-LAT is sensitive to gamma-rays with an energy up to 



approximately 300 GeV, which is too low for the purpose of detection of the monoenergetic 
gamma-rays studied in this work. On the other hand, IACT telescopes reach energies up to 
about 10 TeV, and for these searches the limit is instead set by the finite energy-resolution, 
which may broaden the gamma-ray peak and make it indistinguishable from the continuous 
background. 

The rest of the paper is organized as follows. In Sec. [Til we introduce the UED model, 
including non-minimal versions and implications for the DM phenomenology. Then, in 
Sec. IHfl we describe the analytical calculations for the process Z 1 Z 1 — > 77. Next, in 
Sec. IIV[ we present our results for the cross section and the flux, and compare them to the 
case of B l DM. Finally, in Sec. [V] we summarize our results and state our conclusions. In 
addition, in Appendix |AJ we list the Feynman diagrams for the process, and in Appendix [Bl 
we present a number of Feynman rules that were used in our calculations. 



II. THE UED MODEL 

In the UED model, all of the SM fields are promoted to a higher-dimensional flat space- 
time, and hence, from the four-dimensional viewpoint, they acquire towers of KK particles. 
The model is specified by the geometry of the internal space. We consider a single extra 
dimension with radius R, compactified on the flat orbifold S l /Z 2 , which is the simplest ex- 
ample of a phenomenologically viable UED model The orbifold is constructed from the 
circle S 1 by identifying points that are connected by the orbifold transformation y — > —y. In 
order for the action to respect the Z2 symmetry, the five- dimensional fields have to be either 
even or odd under this transformation, which means that half of the degrees of freedom are 
removed compared to a full circle S 1 . 

The orbifolding breaks translational invariance along the fifth dimension, and hence, the 
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conservation of extra-dimensional momentum. However, a discrete symmetry known as KK 
parity, defined as Pkk = (— l) n , where n G N is the KK number, is conserved. The result 
is that the LKP is stable, in analogy to the lightest supersymmetric particle in models with 
R parity conservation. 

In a five- dimensional spacetime, the four- component Dirac spinors are irreducible, i.e., 
there are no chiral spinors. In order for the model to reproduce the SM phenomenology at 
low-energy scales, the would-be chiral components of the spinors are assigned opposite parity 
under y — » —y, i.e., one is even and the other is odd. Since odd fields have no zero-modes, 
this assignment results in Dirac spinors that are effectively chiral for energies below Br 1 . 
At each non-zero KK level, the chiral fermions in the SM are replaced by Dirac spinors with 
the corresponding quantum numbers. 

Compared to the four- dimensional five-dimensional gauge field Am includes an 

additional component A 5 , which behaves as a scalar in the four- dimensional picture. Re- 
quiring A$ to be odd under y — > —y removes the massless zero-mode field A^ from the 
spectrum. At each excited KK level, A5 is eaten by A^ to generate the KK mass term, in 
complete analogy to the SM Higgs mechanism. However, in non-unitary gauges, these fields 
appear in the loop calculations that are presented in Sec. IIII1 

In general, the orbifold S 1 /Z2 allows for boundary-localized terms (BLTs) at the fixed 
points y = and y = ttB. The effect of such terms is to modify the mass spectrum, as well 
as the coupling constants. In particular, they could affect the identity of the LKP. In the 
minimal UED model, it is assumed that all BLTs vanish at the cut-off scale of the model, 
and are only generated at lower-energy scales through renormalization group running j^]. 
Within this model, the only possible WIMP DM candidates is the first KK excitation of the 
U(1)y gau ge b oson, the B 1 . 



In Ref. 



14j , it was shown, using a restricted set of BLTs, that it is also possible for the 



first KK excitations of the Z boson or the Higgs field to be the LKP. In this work, we are 
interested in the case where the Z 1 is the LKP. In general, the BLTs affect the full mass 
spectrum as well as the coupling constants in a way which is difficult to predict. Here, we 
make the simplifying assumption that all particles at the first excited KK level, except for 
the LKP, have equal masses, given by the mass splitting parameter ^fr\ = Mx/mLKP- Using 
1 However, taking into account the graviton and the full range for the SM Higgs boson mass, the LKP could 
also be the first KK excitation of the graviton or the charged component of the Higgs field 28 1 . 
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similar assumptions, the authors of Ref. 10J performed relic density calculations for the Z 1 
as the LKP, with the conclusion that its mass should be in the range from 1800 GeV to 
2500 GeV. 

For the remaining parameter space, we assume that the coupling constants are not affected 
by the BLTs, which means that the interactions in the non-minimal UED model are identical 
to those in the minimal one. Then, most of the Feynman rules for the gauge self-interactions 
can be found in Ref. 29[. However, a few vertices involving A 5 fields, as well as the vertices 



involving Higgs scalars, had to be derived specifically for this work; those are presented in 
Appendix iBl 

The current experimental bounds on Z 1 DM are quite weak. The strongest direct 
detection limit comes from the XENON100 experiment 30], which constrains the spin- 
independent DM-nucleon cross section. This quantity was calculated for Z 1 DM in Ref. [lo| 
for a set of parameter values that are similar to the ones that are used in this work. Compar- 
ing the results to the XENON100 limits, we find that, for a relative mass splitting between 
the Z 1 and the first-level KK quarks that is larger than a few percent, the model is only 
constrained for Z 1 masses below approximately 1 TeV, which is too small for the correct 
relic abundance. In addition, in Ref. 20[ , it was found that the indirect neutrino signal from 



annihilations of Z 1 DM particles in the Sun is too weak to be observable in current neutrino 
telescopes. 



III. THE PROCESS Z 1 Z 1 77 

As there exist no tree-level Feynman diagrams for the process Z l Z l — > 77, the lowest- 
order contributions are one-loop diagrams. In our calculations, we take into account only 
the contributions from the first excited KK level, as higher levels are suppressed by the 
increasing mass scale. Furthermore, we neglect electroweak symmetry breaking, since the 
mass scale for the KK modes is much larger than the vacuum expectation value of the Higgs 
field. This means that all zero-mode particles are effectively taken to be massless. 

The main differences between the two-photon annihilation processes for B 1 and Z l DM 
is that the Z l , being a non-Abelian gauge boson, has self-interactions with the charged 
SU(2)l gauge bosons, and that it only couples to doublet fermions. The latter point has 
the effect that the contribution from the fermion sector is reduced by a factor of about 
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three. The gauge self-interactions give rise to a large number of relevant Feynman diagrams 



diagrams have 
and used it 



2l|. As a 



that do not exist for the analogous B 1 process. In order to ensure that al 
been included, we have implemented the model in the LanHEP package 
to generate model files for FeynArts [321 ]. which we in turn have applied to generate the 
complete set of diagrams for the process. The full set is presented in Appendix |A 

For the analytical calculations, we closely follow the method described in Ref. 
check of our calculations, we have been able to exactly reproduce the analytical results of 
that work for the B l annihilation process. The amplitude for the process can be written as 

At(Pl,P2,P3,P4) = 4{Pl) e 2(P2)4(Pz) e l(P±) M^pa(pi,P2,P3,P4), (1) 

where the €i{pA are the gauge boson polarization vectors. In the approximation that the DM 
particles annihilate at rest, we have p\ = p 2 = p = (m^p, 0), and hence, 2p + p 3 + p 4 = 0. 
Using this relation, as well as the transversality of the polarization vectors, the amplitude 
can be written in the general form 



M^( Pl ,p 2 , P3 ,p,) 
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where A, Bi, and Cj are dimensionless coefficients that depend on the external momenta 
Pi. Equation ([2]) can be simplified further by using the symmetry for the amplitude of a 
process with identical bosonic external particles, which gives the following relations for the 
coefficients 

B 2 = -B 4 , (3) 
B 3 = -B s , (4) 
B 2 = ~B 3 , (5) 
C 2 = C 3 . (6) 

Indeed, these relations serve as a consistency check of our results. 



The gauge invariance of the model has been studied in detail in Refs. 33|, |34j, where 



it is concluded that the five- dimensional gauge invariance remains in the four-dimensional 
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theory, provided that the KK expansion of the Yang-Mills Lagrangian is performed in the 
correct way. However, the higher-dimensional theory is a non-renormalizable effective theory, 
requiring an ultraviolet (UV) completion at some high-energy scale. Hence, the amplitude 
for the process can generally not be expected to be explicitly finite within the framework of 
the effective theory. The details of the cancellation of divergences depend on the unknown 
UV completion, and cannot be explicitly demonstrated without detailed knowledge of this 
high-energy sector. Thus, we consider only the finite parts of the amplitude, using the fact 
that the total amplitude has to be finite in any renormalizable theory. For a more detailed 
discussion on this point, see Appendix For the same reason, we cannot use the Ward 
identity to further simplify the form of the amplitude, as was performed in Ref. 21]. 

As a partial check of our implementation of the calculations, we have calculated the 
amplitude for the similar SM process ZZ — > 77 in four dimensions, and confirmed that the 
Ward identity and the cancellation of divergences hold for that process. 

Using the Mathematica package LERG-I 35j, we have calculated the coefficients in 
Eq. ()2]) o Due to the large number of contributing Feynman diagrams, the resulting expres- 
sions are complicated and not very illuminating, and hence, we only provide the numerical 
results presented in Sec. IIVI The complete analytical expressions for the contributions from 
the gauge boson self-interactions can be found in Ref. 36j. Finally, the cross section is given 
by 

av = a su(2)«L r 2 + 2 _ 2 + 2 + 2 + 2 
I44irm z zl I 

+2Re [A {B{ + Bl + C*) + B ± B* 6 + 35iC* + 45 6 C* + 2B 6 C; + 6C1C*] } , (7) 

where we have used the bosonic exchange symmetry to reduce the number of independent 
coefficients. 



IV. NUMERICAL RESULTS 



In the numerical calculations, we have employed the same assumptions as in Ref. 



2l| 
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order to enable comparisons to the previously existing results. For the coupling constants, 
2 Note that the software LERG-I 3f| that we have used is suited only for calculations in the Feynman- 
't Hooft gauge, in which the propagators of the gauge bosons have particularly simple forms. However, 
we know from Refs. 33, |3_^ that our setup is indeed gauge invariant, and all the cross checks we have 



performed prove the reliability of our results. 



we have used the values oQ^ = 123 and Ogu( 2 ) = 95, which are the values obtained at 1 
TeV using the ordinary SM renormalization group running of the physical parameters. As 
discussed in Appendix [Cj the accuracy of our predictions depends on the nature of the UV 
completion of the model, which could, in the worst possible case, give a contribution which 
is of the same order of magnitude as our results. Hence, unless the UV contribution would 
almost cancel the contribution that we take into account, our results are valid at least up 
to an 0(1) numerical factor. Note that such a cancellation would require fine-tuning, which 
we have no reason to expect. 

Let us shortly discuss the continuum part of the spectrum, which has been calculated 
in Ref. {3^]. In that work, it has been shown that the continuum spectrum is negligible 
compared to the peak signal investigated here. There are three reasons for this: First, Z 1 Z 1 
pairs mainly annihilate into pairs of W^-bosons {^(J, which contribute insignificantly to the 
continuum spectrum near the peak. Second, the number of diagrams contributing to the 



peak enhances this signal in a significant way, in contrast to the B B case [2l|. Third, 



since the requirement of having the correct DM abundance forces the mass of the Z l LKP 
to be comparatively large, the continuum part of the spectrum is suppressed more strongly 
compared to the peak for larger DM masses {3?]]. Hence, we can safely neglect the continuum 
spectrum in the following, since it is several orders of magnitude smaller than the peak signal 
in the relevant energy range. 

In Fig. [U we show the annihilation rate for the process Z l Z 1 — > 77 as a function of the 
mass splitting between the LKP and the rest of the first-level KK spectrum, ^fr\ = Mi/ttilkp- 
The result is plotted for values of the LKP mass of ttilkp = 1800 GeV, 2250 GeV, and 
2700 GeV, which are the lower end point, the middle point, and the upper end point of the 
range given by the relic density calculations, respectively. As expected, the rate decreases 
with increasing mass, and also with increasing mass splitting. The difference in the rates 
between the highest and the lowest LKP mass is about a factor of two. 

The total annihilation rate as a function of the LKP mass, for the Z 1 as well as the B 1 , is 
shown in Fig. The solid sections of the curves indicate the mass ranges that reproduce the 



relic density for the respective LKP candidates, using the results from Ref. [10j|. Comparing 
the two options in their relevant mass ranges, it is interesting to note that the rate for Z 1 
annihilations is larger than the rate for B 1 by roughly one order of magnitude, even though 
av ~ l/ m LKP anc ^ ^ e mass ran g e f° r Z 1 is considerably higher. The reason is simply 
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FIG. 1. The annihilation rate for the process Z Z — > 77 as a function of the splitting between 
the LKP mass and the rest of the first- level excited KK spectrum, y/rj = Ml/tolkp- 

that, due to the Z l being a non-Abelian gauge boson, the number of contributing diagrams 
increases drastically, therefore leading to a significant enhancement of the annihilation cross 
section. In addition to the 1 / "t^lkp behavior, there is a small mass-scale dependence in the 
Z 1 annihilation rate due to the proportionality of some of the vertex factors to R" 1 (see 
Appendix IB]) . 
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FIG. 2. The annihilation cross-section for KKDM as a function of the mass, for the mass splitting 
Mi/mLKP = 1.1. The solid parts of the curves indicate the mass ranges that give rise to the correct 
relic abundance for the respective DM candidates. 
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Now, the flux at Earth is given by 



*(An) = (2.92 ■ KT 11 nU 2 s^ 1 ) ( ™ 3 ) (M™} ( j GC ) A , An , (8) 

V ID za cm' 3 s 1 / \ mLKP / 

where Af2 is the angular acceptance of the detector and («/gc)aq is the averaged line-of-sight 
integral, which depends only on the DM distribution. We assume a Navarro-Frenk-White 
(NFW) halo profile [39J with a = 1.0, p = 3.0, 7 = 1.0, and r s = 20 kpc. Furthermore, 
we take AQ = 10~ 5 , giving (Jqc) au^Q = 0.136, where b is a boost factor, which may arise 
if the DM distribution is clumpy rather than smooth (^J. In the experimental data, the 
peak will be broadened due to the finite energy resolution of any detector. This effect can 
be represented by a smearing of the differential flux by a Gaussian distribution function. 
Due to the additional suppression by 1/^lkp i n Eq. ([8]), the flux from B l particles is of 
the same order of magnitude as that from Z 1 particles, when they are compared for their 
respective typical masses. Numerically, we obtain the fluxes Qz 1 — 1-2 ■ 10 _12 m _2 s _1 for 
m z i = 2250 GeV and $ B i ~ 1.1 • lO^m^s" 1 for m B i = 800 GeV. 

Finally, in Fig. [3l we show the integrated flux above the threshold energy £ls £i function 
of E~[, and compare the results to existing experimental limits and sensitivities for future 
experiments. Note that, since the contribution from the continuum spectrum is large only 
at low-energy scales, the integrated flux above energies close to m LKP are not affected by 
this contribution. In the mass range that is of interest for this work, the current strongest 
constraints are given by the H.E.S.S. collaboration 23[. Note that the mass range around 1 
TeV is where H.E.S.S. has the best sensitivity. With a boost factor of only b ~ 300, the flux 
for Z 1 DM reaches the limit set by H.E.S.S. In the future, stronger constraints are expected 
to be set by the Cherenkov Telescope Array (CTA) 41]. These limits are reached already 



with a boost factor b ~ 40 for the Z , whereas b ~ 100 would be needed for the B 1 . 



V. SUMMARY AND CONCLUSIONS 



We have investigated the annihilation of Z 1 DM into monoenergetic gamma-rays, and in 
addition, we have confirmed the results for the corresponding process for B 1 DM that already 
exists in the literature. Although the B 1 is the most frequently studied DM candidate in 
the context of KKDM, and the one that is obtained in the minimal UED model, the Z l has 
been shown to be a viable DM candidate in non-minimal UED models. 
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FIG. 3. The predicted integrated flux above the threshold energy E~, with an experimental energy 



resolution given by a, as a function of E~ for the Z 1 and the B 1 , as well as the 



by a number of experiments. The experimental data is taken from Refs. 



42 



43| 



upper limits given 



Our results illustrate that, even when taking into account the expected larger mass of 
the Z 1 compared to the B 1 , the annihilation rate for Z l particles is about one order of 
magnitude larger than the corresponding rate for B 1 particles. However, the flux at Earth, 
which is further suppressed by the square of the mass, is larger for the B , although only 
by a small factor. 

The experimental signature of the process would be a peak in the differential flux on top 
of the continuous gamma-ray spectrum from DM annihilations, which extends up to the 
energy scale ttilkp- The prospects for detection depend on the size of the monoenergetic 



peak in relation to the continuum flux. As was shown in Ref. [2l|, the characteristic peak 
could be effectively erased if it is too low and/or the experimental resolution is not good 
enough. 

It is clear from relic abundance calculations that an experiment reaching up to several 
TeV in energy would be needed to detect the peak. We have shown that the contribution 
from only the monoenergetic flux would reach the existing H.E.S.S. limits or upcoming CTA 
limits already with a moderately large boost factor for Z l DM, whereas a somewhat larger 
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boost factor would be needed for B l DM. 
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Appendix A: Feynman diagrams 

We have identified a total number of 88 Feynman diagrams that contribute to the process 
Z Y Z X — > 77 at one-loop level. In Figs. HI and [61 we show the diagrams with internal 
fermions, gauge fields, and Higgs scalars, respectively. 




FIG. 4. The one-loop level Feynman diagrams involving internal fermion fields for the process 
Z 1 Z 1 — > 77. For each diagram, there is also a corresponding one with crossed final states, making 
the total number equal to 6. The notation A^- 0,1 ^ simply indicates that there are two different 
choices for the flow of KK number in the diagram, and the first or second choice should be used 
throughout the diagram. 



Appendix B: Feynman rules 

In this appendix, we present a number of Feynman rules for the UED model that are, to 
our knowledge, previously not existing in the literature. The fields W§ denote linear com- 
binations of the charged SU(2)l gauge bosons that are analogous to their four-dimensional 
counterparts. The ghost fields are defined as = (c 1 =F c 2 )/\/2, i.e., they are not each 
other's Hermitian conjugates. The constant g is the SU(2) L gauge coupling constant, e is 
the electromagnetic coupling constant, and M\ = is the mass scale of the first excited 
KK level. For the Higgs doublet, we use the convention 



*= * , \ ( B1 ) 
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Gauge boson self- interact ions 
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Gauge boson-Higgs scalar interactions 
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Appendix C: Cancellation of divergences 

In this appendix, we argue why potential divergences in our calculation can be ignored. 
In general, using dimensional regularization, a one-loop amplitude can be written in the 
form 



Mi. 



loop 



A; 



(Cl) 



where Aj/e are divergences and Fi are the finite parts. Per definition, a theory is (pertur- 
batively) renormalizable at the one-loop level if all divergences appearing in the one-loop 
diagrams can be absorbed in counter terms, which must exist as explicit tree-level vertices 
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in the Lagrangian. Hence, for any process that lacks such a tree-level vertex (such as the 
annihilation Z X Z X — > 77 that is investigated in this paper) the coefficients Aj must sum to 
zero. 

Now, if we consider only a finite subset of diagrams, e.g., taking into account only the 
contributions from the lowest KK-modes, we cannot expect that the divergent parts will 
cancel. Although, within this approximation, the interactions are given by a finite number of 
operators which would seem to be renormalizable by dimensional analysis, a complete gauge 
theory is generally not obtained. This means that the cancellation mechanisms arising from 
gauge invariance, i.e., the Ward identity, do not work out. On the other hand, if we sum 
all diagrams at a certain order in the coupling constant, for the complete infinite KK-tower, 
we will recover the full five-dimensional theory at this order in the coupling. However, the 
full five-dimensional theory is non- renormalizable sj], since all coupling constants in the 
five- dimensional Lagrangian have negative mass dimensions. Hence, there is no reason to 
expect that J2i Aj = in this case either. 

However, this issue is only a reflection of the fact that the UED model is just an effective 
low-energy theory, which must be replaced by a UV-complete theory at some high-energy 
scale. The additional contributions to the amplitude in Eq. fIClj) from the UV part of the 
theory lead to the full 1-loop amplitude 

-Mi-ioop = AW + ^ + F uv . (C2) 
If the UV-complete theory is renormalizable, then it must hold that 

^A, + A uv = 0. (C3) 

i 

Without detailed knowledge of the UV completion, it is not possible to show this cancellation 
explicitly, even in principle. However, consistency of the extra-dimensional framework at 
high-energy scales requires that it does exist and that it is renormalizable, which means that 
the divergences Aj arising from the low-energy part do not spoil the predictions. Although 



there exist examples of such renormalizable UV-completions [44|], the exact form of this 
underlying theory is irrelevant for our purposes. 

The size of the contribution from the UV sector, Fyy, can be estimated to be of the 
order of Mkk/Auv = -R~ 1 /Auv ) where Auv is the energy scale of the UV sector, i.e., the 
cutoff scale of the effective extra-dimensional model. In the five- dimensional UED model, 
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the cutoff scale is of the order of 20R 1 



451 ] . i.e., R 1 /Auv ~ 5 • 10 2 . In the worst-case 



scenario, when the UV sector is strongly coupled, this number could be multiplied by a 
factor which is as large as 4w, in which case -Puv/^i-ioop = ^(1)- 

As for the contributions from higher KK-modes, there are diagrams that conserve KK 
number and diagrams that break conservation of KK number, but preserve KK parity. The 
vertices that break KK number are already loop-suppressed, and hence, the corresponding 
diagrams are expected to give only a small contribution. The diagrams that do conserve 
KK number necessarily contain a larger number of massive propagators than the diagrams 
that we take into account (all of which contain massless propagators), and therefore, these 
diagrams are also suppressed, due to the minimum size of 1/R. Quantitatively, the effect of 
higher KK-modes was numerically investigated in Ref. 2l| for the case of internal fermions 
in the B 1 B 1 — > 77 process. In that work, it was found that the contribution is only a few 
percent, and we expect a similar result for our case. 

To summarize the effects of the finite contributions that we neglect, the contribution from 
the higher KK-modes is negligible, while the UV sector could possibly give a contribution as 
large as 0(1), if it is strongly coupled. In that case, our results are valid only up to an order 
of magnitude. On the other hand, if the UV sector is weakly coupled, its contributions are 
suppressed by powers of small coupling constants, and then, our results are more accurate. 

Finally, any loop-calculation involving gauge theories in UEDs where a cancellation of 
divergences does take place at each KK level has the feature to resemble some renormalizable 
theory. For the process B^B 1 — > 77, the Ward identity and the cancellation of divergences 
do hold at each KK level, as was confirmed in Ref. 21|. This is a natural consequence of the 



fact that this process, at any KK level, is essentially a copy of the corresponding process in 
ordinary Quantum Electrodynamics (QED). Hence, the cancellation of divergences in QED 
simply translates to the UED process level-by-level. The same statement is true for the 
fermionic part of the process Z 1 Z 1 — > 77, which our calculations confirm, but not for the 
full process. 
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FIG. 5. The one- loop level Feynman diagrams involving internal gauge fields for the process 
Z 1 ^ 1 — > 77. For each diagram in the first four rows, there is also a corresponding one with crossed 
final states, making the total number equal to 60. The notation A^ ' 1 ^ simply indicates that there 
are two different choices for the flow of KK number in the diagram, and the first or second choice 
should be used throughout the diagram. 
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FIG. 6. The one-loop level Feynman diagrams involving internal Higgs scalars for the process 
Z X Z X — > 77. Here, x denotes the charged component of the Higgs doublet. For each diagram in 
the first row, there is also a corresponding one with crossed final states, making the total number 
equal to 22. The notation A( 0,1 ) simply indicates that there are two different choices for the flow of 
KK number in the diagram, and the first or second choice should be used throughout the diagram. 
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